Astronomy Letters 2009, Vol. 35, No. 10, pp. 688-695. Translated from Pis'ma v Astronomicheskij Zhurnal. 



Effect of the reverse shock on the parameters of the observed 
X-Ray emission during the 1998 outburst of CI Cam 

E. V. Filippova^ '^*, M. G. Revnivtsev^'^, and A. A. Lutovinov^ 
^ Space Research Institute, ul. Profsoyuznaya 84/32, Moscow, 117997 Russia 

^ Max-Planck Institut fiir Astrophysik, Karl Schwarzschild Strasse 1, 86740 Garching-bei-Miinchen, Germany 
^ Excellence Cluster Universe, Technische Universitat Miinchen, James-Franck-Strasse 1, D-85748 
Garching, Germany 



Received: 10 February 2009. 



Abstract. Based on the model of interaction between spherically symmetrical expanding matter and the external 
medium, we have estimated the parameters of the matter heated by the shock that was produced in the envelope 
ejected by the explosion of a classical nova during its interaction with the stellar wind from the optical companion. 
Using this model, we have shown that the matter ejected during the outburst in the system CI Cam had no steep 
velocity gradients and that the reverse shock could heat the ejected matter only to a temperature of ^ 0.1 keV. 
Therefore, this matter did not contribute to the mean temperature and luminosity of the system observed in the 
energy range 3-20 keV. 
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1. Introduction 

Classical nova outbursts can be accompanied by emission 
in both the standard (1 - 10 keV) and hard (> 20 keV) 
X-ray energy bands. Matter heated by the shock that was 
produced by a high-velocity (~ 1000 — 4000 km/s) ex- 
pansion of the envelope of a white dwarf is believed to 
be the source of this X-ray emission. In the system CI 
Camelopardalis (CI Cam), the optical companion is a B[c]- 
type B4III-V star (Barsukova et al. 2006) with a strong 
stellar wind (Robinson et al. 2002; Filippova ct al. 2008), 
which produces a dense circumstellar medium around the 
white dwarf. Filippova et al. (2008) (hereafter Paper I) 
showed that, in this shock is generated by the enve- 

lope expansion in the circumstellar medium, which could 
heat a large amount of stellar wind matter up to 10-20 
keV sufficient to produce a high X-ray luminosity. 

In the envelope itself, a shock (or initially a reverse 
rarefaction wave that will transform into a reverse shock 
as the envelope expands) will also form. Under certain 
conditions, this shock can also heat the matter to high 
temperatures. For example, it follows from the analytical 
calculations by Chevalier (1982) and Nadyozhin (1985) 
that when the ejected matter interacts with a constant- 
density medium, the temperature at the reverse shock for 
an power in the ejected matter density profile of ~ 6 — 8 
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can be lower than the temperature at the forward shock 
only by a factor of ~ 2 — 4. 

A schematic view of the system of shocks produced 
by the interaction of ejected matter with the circumstel- 
lar medium is shown in Fig. 1 (the outer boundary of the 
envelope is a contact discontinuity). In general, the for- 
mation time and the law of motion of the reverse shock 
depend on the matter velocity and pressure distributions 
inside the envelope. However, at present, there is no com- 
plete model that would consistently describe the evolution 
of the profiles of these parameters in the matter ejected 
by the explosion of a classical nova (see, e.g., a review of 
theoretical models in Fricdjung (2002)). On the basis of 
numerical calculations and their comparison with obser- 
vational data during a classical nova explosion, two mech- 
anisms of matter ejection at the initial time are suggested: 
under the action of thermal pressure and through a shock 
wave. The consequences of the ejection of matter through 
these mechanisms were considered by Sparks (1969). He 
showed that in the case of pressure-driven expansion, the 
envelope has a very shallow matter velocity gradient, while 
in the case of shock-driven ejection of matter, it expands 
with a steep velocity gradient. At later expansion stages, 
after the maximum optical brightness of the nova, the ex- 
pansion of matter is described by the model of an opti- 
cally thick wind (Kato and Hachisu 1994; Hauschildt et 
al. 1994). There also exist theories predicting that shortly 
after the maximum optical brightness, the velocities of the 
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envelope layers closer to the white dwarf are higher than 
those of the outer layers (McLaughlin et al. 1947, 1964). 

An example of the possible development of a classi- 
cal nova explosion was given by Prialnik (1986), who cal- 
culated a complete cycle of the evolution of a classical 
nova explosion, from the phase of accretion to its resump- 
tion; the conditions for the generation of a shock wave 
in the envelope were met. According to these calculations, 
within the first half an hour (~ 2000 s) after the thermonu- 
clear explosion, the white dwarf photosphere expands to 
~ 10i?o due to the shock breakout. Within the next ~ 4 h, 
the envelope expands even more due to the radiation pres- 
sure from the white dwarf surface, with the optical flux 
from the system reaching its maximum. It follows from 
observations that the characteristic time it takes for the 
optical flux to reach its maximum for most classical novae 
is < 3 days, but exceptions are also observed; for exam- 
ple, during the outburst of Nova LMC 1991 this time was 
> 13 days (Schwarz et al. (2001) and references therein). 
According to calculations, the outer layers of the envelope 
at this time expand with constant velocities, which in- 
crease toward the outer boundary and reach ~ 3800 km/s 
at it. This part of the envelope ceases to be connected 
with the white dwarf and expands by inertia, interacting 
only with the circumstellar medium. The remaining part 
of the envelope initially contracts under the gravitational 
force and, after some time, again begins to expand under 
radiation pressure in the regime of an optically thick wind. 

It was shown in Paper I that the main peculiarities of 
the behavior of the light curve and radiation temperature 
during the X-ray outburst of CI Cam in the first (spher- 
ically symmetric) approximation could be described in 
terms of the radiation model of stellar wind matter heated 
by the forward shock produced in a classical nova explo- 
sion. In this model, the envelope is ejected from the white 
dwarf due to explosive thermonuclear burning, which al- 
ready on 0.1-0.5 day after the explosion onset has an 
expansion velocity of ^ 2700 km/s and flies under the 
action of an external force, for example, the radiation 
pressure from the white dwarf, with a constant velocity 
for the first ^ 1 — 1.5 days. Subsequently, the envelope 
probably becomes transparent and decelerates, interact- 
ing with the matter of the stellar wind from the optical 
companion. Based on a comparison of the observed rise 
in luminosity with the theoretical dependence, we esti- 
mated the stellar wind density near the white dwarf to 

be no{r < rc) ^ 8.6 x 10^(i2fepct^27oo^'^'^~'^j which trans- 
forms into the law uq ^ at r > Tc = 1.9 x 10^"^ cm. 
In the simplest model, this stellar wind density distribu- 
tion corresponds to a mass loss rate of the optical star 
^ (1 — 2) X 10~^ MQ/yr. The observed time dependence 
of the temperature of the emitting matter at late enve- 
lope expansion stages allowed us to constrain the mass of 
the ejected envelope based on our model, 10"'' — 10~^Mo. 
Note that in this model, the processes in the envelope itself 
were disregarded; in our calculations, we used a finite-mass 
piston as the envelope. 



In this paper, we calculated the contribution from the 
emission of the ejected envelope matter heated by the re- 
verse shock to the observed radiation temperature and 
luminosity of CI Cam during its X-ray outburst in 1998. 




Forward Shockwave 



Fig. 1. Scheme of the shocks produced by the interaction 
of the expanding envelope with the circumstellar medium: 
U is the velocity of the contact discontinuity or the enve- 
lope (depending on the model considered) and D is the 
velocity of the forward shock. 



2. Numerical calculations 

To model the reverse shock, we used a numerical scheme 
described in Paper I: a one-dimensional, spherically sym- 
metric code in Lagrangian coordinates with a staggered 
mesh (the cell radius, velocity, and mass are determined 
at the cell boundaries, while the density, pressure, and 
internal energy are determined at the cell centers). 

At the initial time, the outer boundary of the envelope 
was placed at a distance of 10^^ cm. The density of the 
circumstellar medium was specified as follows: no = 8 x 
lO^CTO""^ at r < Tc and no ~ at r > rc- The initial 
cell size was Ar = 10^° cm. The velocity of the matter at 
the inner and outer boundaries was specified by a time- 
independent constant. 

As in Paper I, we took into account the radiative cool- 
ing of the matter heated by the shocks in an optically thin 
regime. 

As was shown in Paper I, the matter behind the 
shocks is a multitemperaturc plasma in the sense that 
the plasma temperature is nonuniform along the radius. 
Consequently, the radiation temperature that we measure 
based on X-ray observations is an average quantity and it 
may not be equal to the temperature at the shock front. 
Therefore, to obtain the calculated mean temperature, we 
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used the same averaging procedure as that for observa- 
tions. We calculated the ratio of the fluxes in the 3-5 and 
5-20 keV energy bands, which, in turn, corresponds to a 
certain temperature in the radiation model of a single tem- 
perature, optically thin plasma. The method is described 
in more detail in Paper I. 

2.1. Effect of the matter temperature in the envelope 
on the calculated parameters 

The temperature or pressure of the envelope matter is a 
parameter that, in general, can affect the formation and 
propagation of shocks after the decay of an arbitrary dis- 
continuity. 

To understand what the temperature and pressure dis- 
tributions in the ejected envelope are, we can turn to ac- 
tual observations of classical novae. Two outbursts of no- 
vae (Cyg 1992 and LMC 1991) that were observed be- 
fore the maximum optical brightness was reached and 
for which the radiation temperature was measured are 
known to date. However, the effective radiation tempera- 
ture obtained in such an analysis of observations is not a 
good indicator of the physical temperature in the envelope 
(Hauschildt et al. 1994; Schwarz et al. 2001). Nevertheless, 
since the emission from the envelope matter has a max- 
imum in the ultraviolet, we may assert that the matter 
temperature does not exceed ~ 0.1 keV in order of mag- 
nitude. 

To answer the question of how the envelope matter 
temperature affects the propagation of shock waves, we 
performed calculations with the following initial param- 
eters of the matter in the envelope: the density is con- 
stant along the radius (the envelope mass was lO~^A/0, 
the velocity is also constant along the radius and equal 
to 2700 km/s, and we considered several matter tempera- 
tures: 10^,10'', and 10^ K. 

Figure 2 (left panel) shows the matter temperature 
profiles in the envelope and stellar wind on 0.8 day af- 
ter the onset of expansion: the dotted, dashed, and solid 
lines correspond to T = 10^, 10*, and 10'^ K, respec- 
tively. It follows from the figure that the expected range 
of envelope matter temperatures affects weakly the prop- 
agation dynamics and strength of the forward shock and 
leads to unimportant differences in the time dependence 
of the mean radiation temperature (Fig. 2 (right panel)). 
Therefore, below, the initial temperature of the envelope 
matter in our calculations was set equal to a constant, 10'* 
K, unless stated otherwise. 

3. Transformation of the reverse rarefaction wave 
into a reverse shock 

As was noted in Paper I, the interaction of the ejected 
envelope with the circumstellar medium can give rise to a 
reverse rarefaction wave at the very outset. The instant of 
the subsequent transformation of the reverse rarefaction 
wave into a reverse shock depends on the radial distribu- 
tion of envelope matter parameters (such as the velocity 



and pressure). In the same paper, we made very simple 
estimates of the conditions under which the reverse rar- 
efaction wave is generated and the time when it trans- 
forms into a reverse shock. It follows from these estimates 
that the reverse rarefaction wave transforms into a shock 
almost immediately. 

For a clear demonstration of this phenomenon, we per- 
formed calculations in which this transformation could be 
traced in more detail. It should be noted that the param- 
eters specified as the initial conditions bear no relation to 
the actual values: for example, in order that the reverse 
rarefaction wave could recede noticeably from the contact 
discontinuity, the matter temperature at the outer bound- 
ary of the envelope was set equal to 10 keV, but in order 
that the disturbances arising at the inner boundary have 
no time during the calculations to propagate over the en- 
tire envelope, the matter temperature was specified by a 
linear function of the radius and was ^ 0.01 keV at the 
inner boundary; the energy losses through radiation were 
disregarded. The expansion velocity of the envelope mat- 
ter was 600 km/s. 

Figure 3 shows the temperature and pressure profiles 
obtained in this model in the interacting region at various 
times. We clearly see how the rarefaction wave is gener- 
ated (solid line) and how it transforms into a reverse shock 
(in the profiles drawn by the dashed line, the reverse shock 
is seen clearly). The cell number is along the X axis, with 
the contact boundary being located on cell no. 800; the 
forward shock propagates through the stellar wind right- 
ward; the rarefaction wave propagates through the ejected 
matter leftward and transforms into a shock. 

4. Homologous envelope expansion ~ r 

Classical nova explosions resemble in mechanism type la 
supernova (SN la) explosions, thermonuclear explosions 
of white dwarfs (Woosley and Weaver 1986). However, in 
the former case, the explosion energy and, hence, the ki- 
netic energy of the ejected envelope are much lower than 
those in the latter case. For classical nova explosions and 
supcrnovae, the kinetic energy of the ejected envelope is 
estimated to be - 10'*'^ - lO''^ erg (Starrfield et al. 1976) 
and ^ 10^* erg (Khokhlov ct al. 1993), respectively. 

Numerical calculations show that a homologous expan- 
sion of the ejected matter during an SN la explosion is es- 
tablished in ^ 10 s (Ropke et al. (2005) and references 
therein). Dwarkadas et al. (1998) provided the density 
profiles of the ejected matter for several SN la explosion 
models used to describe the observational data. They also 
showed that the density profiles could be described both 
by a power law with an exponent of 7 (nevertheless, the 
exponent can often be different from 7) and by an expo- 
nential law p ^ e~"/"°. 

Chevaher (1982) and Nadyozhin (1985) provided self- 
similar solutions for the decay of an arbitrary disconti- 
nuity in the case of a spherically symmetric homologous 
envelope expansion with powers in the density distribu- 
tion p > 5 into a constant-density external medium. It 
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Fig. 2. Left panel: temperature profiles in the envelope and stellar wind on 0.8 day after the explosion for various 
initial envelope matter temperatures; the solid, dashed, and dotted lines correspond to 10^, 10^ (in both cases, the 
contact discontinuity is at a distance of ~ 1.85 x 10^"^ cm), and 10® K (the contact discontinuity is at a distance of 
~ 1.86 X 10^^ cm), respectively; right panel: time dependence of the mean temperature for these cases. 




Fig. 3. Left panel: temperature profiles in the ejected matter and stellar wind at various times; the cell number is 
along the X axis, the contact boundary is located on cell no. 800. Right panel: pressure profiles at the same times. 



follows from these solutions that when the reverse shock 
is produced, the temperature at the forward shock de- 
creases with time. For example, for p = 5.4, the time 
dependence of the temperature at the forward shock is 
Ti - (V^mm)"^"/^ (Nadyozhin 1985). For an envelope 
with a mass of 10~^Mq, its kinetic energy is 10^^ erg and 
the density of the matter of the external medium is 8 x 10^ 
cm~^, the time tmin is ^ 9000 s; we then find from the 
formula that in the first day of expansion, the temperature 
at the forward shock will fall by a factor of ~ 10, which 
will lead to a decrease in the mean radiation temperature. 

To understand how the homologous envelope expan- 
sion in our problem affects the behavior of the mean radi- 
ation temperature, we performed the following numerical 
calculations. First, we considered the case where the mass 
loss rate from the white dwarf as a result of the shock 
breakout was constant. The power in the radial density 
distribution of the ejected matter is then p = 3. Since the 
velocity of the ejected matter decreases with decreasing 
radius, only the high-velocity outer layers of the envelope 
are actually of interest for the generation of a forward 
shock and energy estimations for the Sedov phase. For our 



subsequent estimations of the kinetic energy and mass of 
the envelope, we took a lower velocity limit of 2000 km/s. 

We performed our calculations for two models: the 
mass of the matter ejected with a velocity > 2000 km/s 
is ~ 10^^ Mq and ~ 2 x lO^^Mg. The kinetic energy of 
the envelope is Ekin,v>200Q ~ 7 x 10^'^ erg in the former 
case and Ekin,v>20oo ~ 1-2 x lO''^ erg in the latter case. 
In calculating the mass, we used the cosmic abundance of 
the matter. However, since the heavy-element abundance 
in the ejected envelope is believed to be higher than the 
cosmic one (Starrfield et al. 1976; Yaron et al. 2005), the 
obtained mass and energy are lower limits. 

Figure 4 shows the time dependence of the mean tem- 
perature of the matter behind the forward shock for these 
two cases (the solid and dashed lines correspond to masses 
of the matter with a velocity > 2000 km/s of 2 x 1O~^M0 
and 10~^Mq, respectively). In Paper I, we estimated the 
explosion energy required to obtain the observed time de- 
pendence of the mean temperature of the matter at late 
expansion stages, when the shock enters Sedov regime (the 
regime in which the shock "forgets" the details of its for- 
mation and evolves self-similar ly). According to these cal- 
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culations, the energy should be ~ (5-10) x lO''^ erg. Thus, 
a further increase in the mass of the ejected matter in the 
model with a homologous expansion will lead to a higher 
(compared to the observations) radiation temperature on 
4-10 days of expansion. 

In Hauschildt et al. (1994) and Schwarz et al. (2001), 
who investigated the structure of the envelopes ejected by 
the explosions of the classical novae Cyg 1992 and LMC 
1991 in the initial expansion period, the constructed mod- 
els that described best the data, had power-law density 
profiles with powers of 15 and 7, respectively. In our case, 
a steeper density distribution in the outer parts of the 
ejected envelope only compounds the situation, because 
the problem has a constraint on the energy of the outer 
layers of the envelope. The velocity of the matter at the 
outer boundary should remain constant, because it deter- 
mines the temperature at the forward shock. Therefore, 
when the parameter p is varied, the mass of the outer lay- 
ers with a velocity > 2000 km/s should be retained. When 
p increases, this will lead to a decrease in the density at 
the outer boundary of the ejected matter at the initial 
time and, hence, the reverse shock will develop more in- 
tensively for some time. The dotted line in Fig. 4 indicates 
the time dependence of the mean temperature for the case 
where the mass of the matter with a velocity > 2000 km/s 
is 2 X IO-^Mq and p = 15. 

We see from the figure and the above estimates that 
the case with a homologous expansion of the envelope mat- 
ter in our one-dimensional model is in conflict with the 
observational data. 

The absence of a homologous expansion phase during 
the outburst in the system CI Cam can be explained in 
several ways. The conditions for the generation of a shock 
wave in the envelope may have not been met during the 
explosion; therefore, the matter was ejected by thermal 
pressure without any velocity gradient (Sparks 1969). Or 
it is possible that the matter ejected as a result of the 
shock breakout was almost immediately stopped due to 
the high density of the external medium. The shock pro- 
duced by it was damped out almost immediately and the 
observed shock was generated by the subsequent ejection 
of matter due to the radiation pressure. 

5. Envelope expansion with a constant velocity of 
matter 

We performed calculations for three density proflles of the 
matter in the envelope at the initial time: p = const, p ^ 
r~^, and p ^ r~^. The envelope mass was taken to be 
~ IO-^Mq. 

To keep the velocity of the forward shock during the 
decay of an arbitrary discontinuity the same as that when 
the piston envelope is pushed, the velocity of the matter 
in the envelope should be higher than the piston velocity. 
In these calculations, we set it equal to 3000 km/s. 

The derived temperature profiles of the matter behind 
the forward and reverse shocks are shown in Fig. 5. It 
follows from this figure that in all cases the temperature 
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Fig. 4. Time dependence of the mean temperature of the 
matter behind the forward shock for a homologous expan- 
sion of the matter in the envelope: the solid line correspond 
to the model with a mass of the matter in the envelope 
having a velocity > 2000 km/s of 2 x 1O~^M0 and p = 
3; the dashed line corresponds to lO^^M©, the exponent 
in the density distribution of the ejected matter is p = 3; 
the dotted line corresponds to 2 x 10~^Mq, p = 15. 

behind the reverse shock does not exceed 0.1 keV (recall 
that radiative cooling in our model switches on at temper- 
atures > 0.1 keV), i.e., this matter does not contribute to 
the observed fiux in the 3-20 keV energy band and to the 
temperature averaged over the X-ray flux during approx- 
imately the flrst 0.7 day of expansion. 

Investigation of the behavior of the reverse shock at 
later times is hindered by the absence of reliable theoret- 
ical models for the distribution of physical parameters in 
the expanding matter. 

Note that the envelope expanded freely in these cal- 
culations. However, as was shown in Paper I, on the first 
day of expansion, an external force keeping its velocity 
constant should act on the matter in the envelope. It is 
quite possible that the reverse shock will be suppressed 
even more in this case. 

We sec from the figure that the radiative cooling of the 
matter behind the reverse shock is important only for the 
model with a density profile p (dashed line) - the 

temperature profile exhibits a "shelf behind the reverse 
shock. Let us show that the radiative cooling in this case 
also takes place in an optically thin regime. 

5.1. Radiative cooling of tlie matter behind ttie reverse 
shocl< 

To establish the regime of radiation of the matter behind 
the reverse shock, we calculated its optical depth. The op- 
tical depth for Thomson scattering up to the region behind 
the reverse shock in which the radiative cooling becomes 
important is less than unity. A significant contribution to 
the absorption of radiation in the envelope matter can be 
made by absorption in lines. However, using the opacity 
tables calculated with the OPAL code, we found that for 
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Fig. 5. Temperature profiles behind the forward and re- 
verse shocks on 0.5 - 0.7 day after the explosion onset for 
an envelope with a mass of - IQ-^Mq: the solid, dotted, 
and dashed lines correspond to p = const, p r^"^ , and 
p ^ r~'^, respectively. The contact boundary is on cell no. 
100. 



a density of the order of 3 x 10^^ cm ^ and a temperature 
of 5 X lO'^ — 10'' K, the absorption cross section in the 



matter doesn't exceed 



10" 



cm . The optical depth 



corresponding to this cross section is also less than unity. 
Thus, the matter behind the reverse shock, just as behind 
the forward shock, radiates in an optically thin regime. 

As we have already said above, the heavy-element 
abundance in the envelope matter ejected by the explo- 
sion of a classical nova can be higher than the solar one, 
which, in turn, can lead to an increase in the cooling rate of 
the matter. In our calculations, we retained the solar abun- 
dance, because the characteristic cooling time for the den- 
sity obtained in the calculations (3x 10^^ cm~'^) is Trad ^ 4 
s even for the solar abundance. If it is lower by a factor of 
several, then this will not affect the calculations. It follows 
from the formulas in Paper I that for the matter behind 
the reverse shock, the characteristic time it takes for a 
Maxwellian velocity distribution of ions at Ti = 0.1 keV 
to be established is tu ^ 6 x lO""*^ s (in these estimates, we 
took the parameters Ai = 1, ^ 1/1836, InA = 15, Zi = 
1, and Ui = rie, as in Paper I). Formally, for the ion and 
electron temperatures Ti = 0.1 keV and Te = 1 eV, the 
ion Maxwellization time is longer than the time of tem- 
perature equalization between the ions and electrons by a 
factor of Tu/Tie = 4.5 x 10-^(Te/Ti + 5x IQ-^)-^/^ ~ 42. 
This means that the ions initially transfer their energy to 
the electrons without still having a Maxwellian velocity 
distribution, while Spitzer's formula for Tie begins to work 
only when > QA3Ti 0.01 kcV. For our estimations, 
we take = 0.05 keV (in this case. Tee ~ 5 x 10~^ s), then 
Tie ~ 5 X 10"'^ s. Clearly, the time of electron heating by 
ions to Te = 0.05 keV should be of the order of the value 
of Tie obtained. The time it takes for an ionization equi- 
librium to be established is Teq ~ Is. Thus, it follows from 
our estimates that the characteristic times it take for an 
equilibrium to be established behind the reverse shock are 



( |http://hea-www.harvard.edu/APEC/REF[ ) to calculate 
the plasma energy loss rate are metQ. 



6. Conclusions 

We investigated the effect of the reverse shock on the ob- 
served parameters of the X-ray emission during the 1998 
outburst of CI Cam using a spherically symmetric model 
for the interaction of the envelope ejected by the nova ex- 
plosion with the circumstcUar matter. Comparison of our 
numerical calculations and observations in the frame of 
this model led us to the following conclusions. 

— Te homologous expansion phase of the matter during 
the explosion in CI Cam most likely was either absent 
or short and did not give rise to an observable forward 
shock in the stellar wind. The velocity profile in the 
matter ejected by the explosion had no steep gradients. 

— For a free envelope expansion at a constant veloc- 
ity and the explosion parameters that we obtained in 
Paper I, the reverse shock could not heat the matter to 
temperatures above ~ 0.1 kcV during the first ~ 0.7 
day of expansion. 

— During the 1998 outburst of CI Cam, the contribution 
from the matter heated by the reverse shock to the 
observed luminosity in the 3-20 keV energy band and 
to the temperature averaged over the X-ray fiux during 
the first ^0.7 day of envelope expansion was negligible 
compared to the contribution from the matter behind 
the forward shock. 
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